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ABSTRAC T 


The results of the March 1971 wind tunnel tests, second series of tests per- 
formed in the NASA Ames 40 x 80 foot wind tunnel, of the DH 2011 jet-flap rotor 
are presented and analysed. The tests have been focused on multicyclic effects 
and the capability of this rotor to reduce the vibratory loads and stresses in the 
blades. 

The reductions of the vibrations and stresses at tip speed ratio of 0.4 have 
attained 50%. The theory shows further reductions possible, reaching 80%. The 
results show that the performance characteristics after the modifications intro- 
duced since 1965 remained unchanged. The domain of investigation has been 
enlarged to include the tip speed ratios of 0. 6 and 0. 7. 

To analyze the complex aeroelastic phenomena a new analytical technique has 
been utilized to represent the mathematical model of the rotor. This technique, 
based on transfer matrices and transfer functions, appears very simple and it 
is believed that this analysis is applicable to many kinds of investigations in- 
volving large numbers of variables. 
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March 1971 Wind Tunnel Tests 



of the Dorand DH 2011 Jet Flap rotor 



Final Report 



INTRODUCTION 


The experimental studies concerning the DH 2011 jet flap rotor have started 
in 1964 by a series of whirl tests outside tunnel (Ref. 1). The first series of 
wind tunnel tests took place in 1965 in the NASA, Ames 40 - by 80 - foot wind 
tunnel (Ref. 2 and 5). The second series of wind tunnel test which related in 
the present report has been preceded by static tests (Ref. 3 and 6). 

The application of the jet flap to helicopter rotor blades has been analytically 
treated in Ref. 4. The practical aspects of the use of jet flap rotors to rotary 
wing aircraft has been investigated in Ref. 7. 

The results of the tests performed in 1965 in the 40 x 80 foot Ames Wind 
Tunnel have showed the very high lift and high forward speed potential for this 
type of rotor. Advance ratios up to 0.5 were reached without indications of re- 
treating blade stall, and the blade loadings i exceeded conventional rotor capabilities 
by factors of 2 and more (Ref. 5). Upon completion of the 1965 tests, it was deci- 
ded to increase the range of jet deflections and to modify the flap control linkages 
so that further investigations could be made into the multicyclic potential for 
relieving vibratory loads and stresses. In fact the DH 2011 rotor possesses the 
unique feature of introduction of multicyclic control effects by use of mechanical 
cams acting on the jet flap deflectors, and the exploration of this capability 
became the main theme of experimental investigation. 

This document presents also the analytical technique developed specifically 
to examine the wind-tunnel results. Ref. 6. In fact the multicyclic effects involve 
a large number of variables (harmonic components of stresses, loads, control 
signals) which require special approach to obtain a clear picture of the phenome- 
na studied. As shown later, the technique used is based on transfer matrices 
and transfer functions and is believed to be more generally applicable to analyze 
mathematically complex models. 
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The DH 2011, a full scale 40 ft in diameter, 4.5% solidity, two bladed rotor 
has been fully described in Ref. 1 through 6. The modifications that have been 
introduced since the first series of tests in the 40 x 80 Ames Wind Tunnel are 
three : 

1. The pneumatic control system has been replaced by a mechanical one to 
control with more precision the multicyclic jet flap deflections, which are 
imposed by predetermined cams. 

2. The amplitude of jet flap deflections has been increased from + 25° to + 50°. 
This increase lias been motivated by desire of higher efficiency of jet flap 
in translational rotor configurations at high advance ratios. 

3. The nozzle height has been increased from 4.2 mm to 5.3 to obtain some- 
what higher power input and also to be able to judge the influence of the 
nozzle height on the general performance of the rotor. 

The initial test program, Ref. 8, has been performed to a high degree, 
roughly speaking to 80%. Two main reasons Imposed the curtailing of the program, 
limited to three weeks in the wind tunnel : 

1. Lack of power of the air generator, whose pressure ratio was limited to 3. 6 
(most of the tests have been done at 3.4), the highest pressure ratio expected 
being 4. 0 to 4. 2. 

2. Premature wear of the mechanical control system that rendered the test 
procedure more complicated. 

During these tests, the rotor has been running a total of 24 hours and 
30 minutes, representing 299 measurement points distributed over 29 runs. 

An additional run was made to determine the hub drag. 

One of the main interest of these tests has been to make possible the study 
of the effects erf a multicyclic control on a jet flap rotor. Numerous measurement 
points have been made for advance ratios between 0. 25 and 0. 7. Two values of the 
blade pitch angle have been tried : © 0 ? = 5° and 3°. A pronounced influence of 
the multicyclic control has been observed on blade stresses and hub vibratory 
forces, preliminary investigations showing stress reductions of about 50%. 

The test program included originally an optimization of o( s . Unfortunately 
only few points have been obtained due to limitations in the controls amplitude and 
it does not seem that they will be enough for the intended purpose. 
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Aside these tunnel tests, simulations of the multicyclic control have been run 
in cooperation with NASA research team, on the Ames IBM 360 digital computer. 

Dynamic data have been recorded, part by the GIRAVIONS DORAND team, as 
oscillagraph records, part by the Ames team, as oscillograph and magnetic tape 
records. Also static values have been read out from the instruments and written 
down by both teams. These values are presented in Section A in two sets of 
tables. The first set corresponds to direct readings. The second set to quantities 
computed from read values by the Ames program FSAO 230. 

Despite the fact that it had not been possible to obtain all the planned measu- 
rement points the March 71 DH 2011 wind tunnel tests should bring to light many 
important informations about jet flap helicopter rotors and the efficiency of intro- 
duction of a multicyclic deflection law for reducing blade stresses or/and hub 
vibrations. 
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PRIMARY TEST DATA 


1. TEST CONFIGURATION PARAMETER TABLES (Table I) 

In these tables are given the values of the parameters defining the 
conditions of the test at each measurement point. These values come from 
direct instrument read-outs. Each test case, is coded by two numbers, the 
first is the run number, the second the measurement point. For example 
the 7th measurement point of run 4 is denoted by ’’run 4-7". This conven- 
tion will be used throughout the report for its convenience. Symbols used are 
those of Ref. 4. Some are restated below. 

The following quantities are given in the tables; 


Symbol 

Meaning 

Units 


Blade pitch angle 

degree 

Vc 

Phasing angle of the multicyclic cam 

degree 

RPM 

Rotor rpm 

r. p. m. 

T 

Rotor thrust 

lb 


Rotor shaft angle (negative forward) 

degree 


Suspension longitudinal tilt 

degree 

Ao 

"Collective" flap deflection 

% of max. 



deflection 

A\ 

"Lateral cyclic pitch" flap deflection 

1 1 


"Longitudinal cyclic pitch" flap deflection 

h 

c? y 

Multicyclic cam amplitude 

% of max. 

J 


amplitude 

Flap deflection peak-to-peak amplitude 

degree 

Ap 

Pressure difference across the flowmeter 

mbar 

pi 

Pressure at the flowmeter input 

bar 

V 

Tunnel velocity 

Ids 


The flap deflection angle is positive downward. Collective and cyclic 
components or given by : 

A 0 - Cos g, Svru H'" 

Maximum values of the imposed deflections are 

max 

A, a 

may 

B 


i max. 


= 50 1 
= 30 1 
= 30 




It should be noted also that 60% amplitude of multicyclic cam 
(i. e. Azimuth = 60) corresponds to the theoretical cam amplitude 
(see paragraph IV for the cams definition). 


n. ROTOR AERODYNAMIC AND PERFORMANCE COEFFICIENT TABLES 
(Table II) 

Basic parameters have been read out for each test point (see 
Table I) and some others (such as lift and drag, derived from the balance 
dynamometers signals) have been obtained by the Ames team. From 
all these measurements, principal aerodynamic coefficients and per- 
formance coefficients have been computed by the Ames team on the 
IBM 36C using the program named FSA 0230. 

These coefficients are given in Table II. They are : 


Symbol Meaning Units 

^ Advance ratio 

Rotor tilt angle degree 

T Thrust lb 

W Air mass flow lb/s 

GrWP Isentropic gas power hp 

feSHp Equivalent shaft power hp 

CJR/(T Rotor iet coefficient 

CLR/<5" Rotor lift coefficient 


G#R /<5" Rotor propulsive force coefficient 
Cy /CT Lateral force coefficient 
CPE / G Rotor shaft power coefficient 
Cpp /C PE Propulsion efficiency ratio 

REMARKS 

1. Jet rotor efficiency 

The efficiency ^ of the jet rotor is defined as the ratio of the avai- 
lable mechanical power on the rotor (ESHP in Table II) to the power 
delivered by the gas generator (GHP). 

*? = 


ESHP 

GIIP 
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It depends upon the pressure ratio Rp and the rotor r. p. m. . The 
pressure ratio is defined here at the flowmeter, i. e. , if P q is the 
static pressure in the tunnel : 


R = 

P 


/ p o 


During the tests, Rp varies 1.7 to 3.6 and is found between 28% 
and 36%. These low values are due to the duct losses (K^ ^ 3) 
and gas temperature which were high in this experimental rotor. 

This feature made also difficult to match the power supply to the 
rotor, and hence, not enough power was available for tests at high 
advance ratio. 


Variation of ^ with Rp are shown in figure 19 for different values 
of the rotor r. p. m. 

Propulsion efficiency ratio 

This quantity is the ratio *?p of the power available for translation 
to the power delivered by the engine (here the gas compressor), 




XV 

ESHP 


or using reduced power coefficients : 


(CPE) 


(see table II) 


The propulsion efficiency ratio is plotted in figure 20 as a function 
of the advance ratio p and the rotor tilt angle oC . For small 
angles, 

^degree ~ ^ X/L 

The curves show increasing values of K P with ^ and demonstrate 
the interest ctf the jet flap rotor for flying at high forward velocities. 


PERFORMANCE COEFFICIENTS PLOTS 

To summarize the numerical results of Table II, different plottings of th- 
values of the most important coefficients are presented here. 


Three coefficients have been considered for different values of Q 


t* ' <* 5 ’ 


and different cam configurations. 
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They are : 


C LR /<$ (lift coefficient) 

C XR (rotor propulsive force coefficient) 

Cj R ^ (rotor jet-m omentum coefficient) 

- All these coefficients are therefore based on an area of 59.4 sqft 
(i. e. 0, 0488 TrR 2 ) which makes them of theCLR/<F type with 

<3= 0. 0488) 

- Plots of C LR vs and C RR vs are given in figures 

1 to 18. 

Figures 1 to 8 correspond to a blade pitch angle ©07 of 5° , each 
plot is made at a given but different . The curves corresponding 
to Ck> & = ct 0 ' result from a fairing of the data points. The same set of curves 
was used in figures 1 to 4, while in figures 5 to 8, a global translation 
along the x axis was introduced in each plot, in an attempt for accounting 
for the change in V/JlR. 

Figure 9 gathers data points corresponding to two values of ^ . 
However, the points are differentiated according to the value of_ the collective 
pitch A 0 in order to obtain a tentative fairing of the curves = cte. 

The same features are found in fig. 9 to 18, corresponding to $ o .= 3°. 
Remarks 


The scatter found in the results was to be expected because some 
other parameters than those mentioned previously had changed from one 
run to another, like for instance the main controls (4,, A„ 6,) and the 
multicyclic controls. This problem is studied in more detailed in Section 
c, m. 

Nethertheless, the general aspect of the performances is in good 
agreement with previous tests of the rotor. 

The average lift coefficient 

C L ® ^ C_r / ^ ^ 6 ( Cl R /<S~ ) 

has reached values above 1, which is characteristic of jet flap rotors, 
but not a limit. Higher values would have been obtained with more avai- 
lable power on the engine. 


MULTICYCLIC CAMS DEFINITION AND REPRESENTATION 
1. Definition of the cams. 


A multicyclic cam may be defined by the polar equation of its section. 
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or more conveniently by its Fourier expansion. For a given blade 
azimuthal angle W , the flap deflection law mechanically imposed 
by the cam, assuming a perfect transmission, is written as : 

1,66 A y. ( CdS 3 c V+ *L) + 4 S'l 2 %) + 4 coS 3 &+%.)+ B, 

+■ • • C .tc J 

(By convention, t) is positive for flap down) 
where 4^ is the angular position of the cam with respect to a fixed 
reference on the hub, and Ay the amplitude of the control (denoted 
as "Azimuth" in Table I). 

For *f c - 0 and A*/* = 60a', A 2 , 6i ... etc. are the Fourier 
coefficient of the theoretical multicyclic flap deflection. 

Four cams have been used in the tests. They are denoted as 
Cams I, III, IV and V, and the corresponding values of An 
and 6 n are given in Table III. 

In general, for given values of and ^c, , the Fourier coef- 
ficients of the multicyclic law, defined by : 

J c 4 S,a a Hf ♦ Jj cos 3V c$ 3 5irt 3 V + a tc I 

are related to the A n and e„ by : 

jcf R = -1,66 ( A a cos + §n. sin. n f±) 

^ c ^ ^ — 1,66 A oj^ 1 cos a ri 4^ ) 

2. Graphical representation of the cams 

It is convenient to represent graphically a given cam by vectors indi- 
cating, for each harmonic, the positions of the blade where the maxi- 
mum flap deflection occurs when only driven by this harmonic. 

This of course will give two opposite vectors for the 2p harmonic, 
three vectors 120° apart for the 3p, four at 90° for the 4p, and 
so on. 


A change Yc. in the phase setting of the cam will only rotate the 
whole pattern by an angle 4L 

In addition, the length of the representative vectors may be made 
proportional to the corresponding harmonic amplitude. 

On figures 21, 22, 23 and 24 are shown this representation of the 
four cams at l f c = 0. 
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The range of the angles swept by the different harmonics of the four 
cams used during the tests are shown on figures 23, 24 and 25. 

The numerical values of the multicyclic components, , of the jet- 
flap deflection, c5 , of five original cams at their 60% control position 
are given below : 


Cam I 
Cam II 
Cam HI 


Cam IV 


Cam V 


- + 11 cos 2 V 

<$y =+&.4 cos 2 V* - 1. 8 cos 3 V" + 1. 3 sin 2 V +4.5 sin 3 V 

dy = 2.5 cos 2 V - 3. 5 cos 3 V + 2. S cos 4V” + 2.3 sin 2 V 
+ 3 sin 3 V 

=2.5 cos 2 V - 3.5 cos 3V + 2.5 cos 4 4^ + $ sin 2V 
-4.0 sin 3 V 

= 7. 0 cos 3 ^ 


The components of are given in degrees. Only four cams have been 
used during tests. The cam N* H has never been used. 
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ANALYSIS OF THE RESULTS 


I. - STRESS AND VIBRATORY FORCE DATA 


During the March 71 wind tunnel test of the DH 2011 jet flap rotor, 
different rotor configurations have been studied with different mult icye lie 
cams. For each cam, several phase settings were used and for a given 
phase setting, the amplitude of the cam was varied (see Table I and Section 
B, IV). 


The oscillograph records taken then, have been analysed at 
GIRAVIONS DORAND to determine the effect of the multicyclic cams. Peak- 
to-peak values cf the stresses and vibrations have been measured and com- 
pared and these results are displayed in the following tables. 

1. - TABLE IV - STRESS REDUCTION (pages 49 to 52 ) 

This table gives the peak-to-peak values of the flap bending stresses 
at 0.45 and 0.7 R. They are denoted respectively as "ampl. 
max 4.2/1" and "ampl. max 4.3/1". 

The column "CAM. AMPL. %" indicates the percentage of total 
multicyclic cam amplitude which was used in the corresponding test 
point. 

The stress reduction is given in % according to the formula : 

&<5 / 6 _ 100 

where 6“ 1 denotes the value of the stress for zero or small cam 
amplitude, <5 2 to an increased value of the amplitude. (Therefore, a 
positive value of &<5 / <5 corresponds to a reduction). 

Most significant cases exhibiting the stress reducing effect of the 
multicyclic cams are plotted in figure 38. 


The run number corresponding to the value <5* 2 is shown on the 
X axis, and the cam type and phase shift used in this case are 
indicated right below. 
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2. - TABLE V - VIBRATORY FORCES REDUCTION (pages 53 to 57 ) 

Peak-to-peak values of the vibratory forces have been measured from 
the oscillograph records, corresponding to the left vertical arm of 
the balance (load cell F 5. 1). 

Conventions for the reduction are the same as for the stresses and 
a similar display of significant cases is given on figure 39. 

Also on figures 28 to 37 are plotted the peak-to-peak values of the 
vibratory forces as a function of the cam amplitude (figures 28 to 
30 contain in addition results for F 5. 2 and F 5.3 load cells 
giving forces in the right aft and forward vertical arm of the ba- 
lance). 

3. - STATISTICAL RESULTS 

Several multicyclic cam settings have been tested in an almost ran- 
dom sweep in order to cover as large as possible the available 
phase-amplitude range for the multicyclic flap deflection (see Section 
B, IV). 

Results of these tests, summarized in table VI below, show that a 
significant reduction has been obtained in about half of the cases. 
Similar results are obtained for the vibrations. 




Stresses 

at 

Aft load 



r=0. 45R r=0. 7R 

cell 



T _ 

T 

F_ 



0.45 

0. 7 

5. 1 

More than 5% 

reduction 

26 

30 

46 

More than 5% 

augmentation 

28 

24 

13 

Less than 5% 

change 

12 

12 

10 

Total number of cases 

66 

66 

69 


TABLE VI - 

Statistical effects of 

the cams 



To display the effect of the cams more systematically, statistical 
effects of the cam phase shift has been studied in the following 
way : 


For a given cam type with a given phase shift, a serie of test 
points had been made in different flight configuration. In each 
flight configuration, a first test point was made with no or small 
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cam amplitude, and the peak-to-peak stress amplitude measured, 
say <5 For an increased amplitude of the cam, an other test 
point was made, and the stress was thenC g. These turn test points 
constitute a ’’case” for which the variation of the stress is com- 
puted according to Eq. 1. If N cases have been considered for this 
cam, an averaged variation is computed as : 

U N 

< / g- ;> --rr?- T ( a gVOL 

This quantity is plotted as a function of the phase shift for each 
cam on figure 40 for the stresses and figure 41 for the vibrations. 
This plottings show the marked effect of the phase and therefore 
the importance of a correct setting for optimum results. 

The analysis ctf the signals has been obtained for different values of Oy. 
The results show that stress reductions are non necessarily correla- 
ted with reduction of vibratory forces. 

Striking examples of stress and vibration reduction obtained by the 
multicyclic control are displayed on figures 42 to 46. Here are 
compared the stress {or vibration force) time histories for diffe- 
rent multicyclic jet flap deflections .which are also plotted. These 
signals have been copied from the original oscillograph records 
taken during the tests. 

As it appears, up to 40% reduction ctf stresses and 48% reduction 
of vibrations have been observed, which demonstrates clearly that 
the multicyclic control of the rotor is a powerful means for redu- 
cing blade stresses and vibrations. 
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II. - ANALYSIS OF THE MULTICYCLIC EFFECTS 


1. - INTRODUCTION 

The overall analysis of the experimental data described in the 
preceeding paragraph has shown that the multicyclic cams have a 
definite effect upon the stresses and vibrations. Due to simulta- 
neous variations in other parameters, this analysis has princi- 
pally a statistical value, but is not sufficient for relating preci- 
sely the stress variations to specific proport ies of the cams. A 
search for a correlation between stresses and flap deflection has 
been motivated by three important factors : 

a) the need to establish the actual contribution of the cam and the 
influence ctf the other parameters, 

b) the application ctf such results to the definition of the best type 
of cam to be used for a systematic stress reduction. 
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c) the possibility of describing this correlation by a linear trans- 
formation as demonstrated on simulated results (Ref. 3). 

Indeed, not only <^s and © 0 . 7 were varied during the tests, 
but, even for a given flight configuration, the introduction of 
multicyclic cams was accompagnied by some changes in the three 
main controls A o , and . Therefore, the first step is to iso- 
late the multicyclic effects and correlate them with the multicyclic 
components of the cams. 

Then the optimization of the cams can be approached in two different 
ways. One is to simulate as well as possible the rotor dynamics 
on a computer and try different types of cams, as this was done 
in the present test. But in this case, the parameters can be set 
exactly as desired and it is possible to bring out the multicyclic 
effects only. This has been done at the AMES Research Center 
using the Evans - Me Cloud program (Ref. 4). For instance, let 
us consider a serie of ten simulations corresponding to the follo- 
wing conditions. 


* 0.7 = 5 ° 

250 r. p. m. 

C JR = 5.4. 10“ 4 


°<s = - 6 ° 
p =0.25 
Clr/ 6 - = 0. 096 


C x p/«t= 0. 007 


Collective pitch 51.1 degrees 

Cyclic (A,) 0 

Cyclic ( s, ) 14.5 " 

The program evaluates the Fourier coefficient of the thrust coeffi- 
cient, C T , which is more or less representative of the vibratory 
forces at the hub. One can then deduce the time history and the 
peak-to-peak values and obtain such table as below : 
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RUN 

CAM 

^2 

COMPONENTS (degree) 
^2 cf 0 V 0 

Ct x 103 
peak-to-peak 

50- 1 

0 

0 

0 

0 

1. 7 

50- 2 

0 

0 

0 

0 

1.7 

50- 3 

- 5 

0 

0 

0 

1.21 

50- 4 

0 

- 10 

0 

0 

2. 04 

50- 5 

- 5 

- 5 

0 

0 

1.24 

50- 6 

0 

0 

- 5 

0 

2. 18 

50- 7 

0 

0 

0 

- 10 

1. 73 

50- 8 

0 

0 

- 5 

- 5 

2. 02 

50- 9 

- 5 

- 10 

- 5 

0 

1.29 

50-10 

- 20 

0 

0 

- 5 

1. 8G 


TABLE VI - Simulated multicyclic effects on vibrations 


From this test, it can be deduced that the cam of run 50-3 is the 
best, but in fact there is no way to tell whether it can be impro- 
ved, particularly by adding a 3 p harmonic component. Actually 
it exists a cam which reduces Cr down to 0.6.10-3, peak-to-peak 
value, as it will be seen later in paragraph III. 5. 

The other problem, is the accuracy of the simulation itself. Com- 
putations based upon an uniform induced flow field is indeed known 
to give good results as far as static values and performances are 
concerned, but disputable ones when higher harmonics are consi- 
dered. 

Since a realistic mathematical model ctf the rotor would be far 
too complicated for evaluating the multicyclic effects, a global 
description of the rotor using its input - output properties has been 
proposed which turned out to be very powerful for anatysing the 
test data. 
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since all the tests were performed in the steady state conditions. 
The original idea, tested first on simulated data, was to relate the 
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Fourier coefficients of the measured stress to those cf the multi- 
cyclic cam, by a linear transformation, considering only the har- 
monics of rank 2 and 3. But with test measurements this method 
turned out to be unapplicable for the following reasons : 

a) Due to some imperfections of the mechanical linkage, the actual 
flap deflection, as measured by the Cimatran, was considerably 
different from the deflection originated at the cam. It became 
therefore necessary to Fourier analyse the Cimatran signal 
instead of using directly the theoretical cam Fourier coefficients 
in the correlation analysis. Comparison between the Fourier 
coefficients of the multicyclic cams ordered deflections and 
those cf the actual flap deflection are given in Table VII for 

30 runs. 

b) For the same reason, runs which were supposed to be perfor- 
med with the same collective and cyclic were not really so. 

Since these controls contribute also to the stresses harmonic 
content, it was thus necessary to include them in the model, 
with the multicyclic deflections. 

With these modifications, partial results were obtained using few 
runs, and they were encouraging, the dispersion on some coeffi- 
cients being less than 10%. However it was soon realized that the 
analysis should at least include the 4th harmonic because the study 
of the Cimatran signals showed that it was of a similar or even 
greater magnitude than the 3 rd harmonic. At this point, a least 
square regression method has been introduced and a 80% correlation 
obtained on 15 test points. 

It appeared then that the value of the rotor shaft angle °^ s , could, 
also be included in the model in a linear way. With this modification 
30 different test points have been correctly represented by the model. 

Finally, this linear representation makes possible the computation 
of a cam which would reduce the stresses. At first an ’’ideal cam" 
was defined according to Ref. 3, which makes zero the 2nd, 3rd 
and 4th harmonic of the stress. But with this particular definition, 
peak- to- peak values were not always reduced. An other criterion 
has been used in recent analyses which leads to a systematic reduc- 
tion in all 30 cases. Ideal cams corresponding to this new criterion 
are easy to compute and could be used as starting point in a finer 
optimization algorithm. 
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For these correlation studies, a general FORTRAN IV program 
has been established- All the test points analysed were made at the 
same advance ratio ( P- 0.4), but provisions have been made 
in this program to include the effect of /* in order to cover a very 
wide range of test conditions. 

PRELIMINARY INVESTIGATIONS 

2. 1. Basic ideas 

In Ref. 3 a method is described for correlating the harmonics 
2 and 3 of the jet flap deflection to those of the thrust coef- 
ficient Crp, or of the moment coefficient M-p. The hypothesis 
made there is, that for all runs used in the correlation ana- 
lysis : 

a) The rotor configuration parameters are the same (i. e. -A 
^JR» ®o.7 an< ^ °^S 

b) The rotor controls are the same (i. e. collective and cyclic 
flap deflection), 

c) The flap deflection does not contain harmonics of rank 
greater than 3. 

The model proposed in this method for describing the rotor 
may be derived simply from the expression of the local 
aerodynamic forces on the blade. Within the assumption of 
small angles (local angle of attack, flapping angle and iet 
flap deflection angle) the force by unit area at the distance 
r of the shaft may be expressed as : 


~ ^ f ( rJi + Vs\n ^ ( 04V O + ( C j + >/cj ) cT 3 


where 


V is the free stream velocity 

Cj the jet momentum coefficient 

S L the jet flap efficiency 

6 the local blade pitch angle 

-f is an aerodynamic angle of attack given by 


rfb ) / ( r -S1_ V <,\ci Y ) 
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with 


V't' induced velocity 

V z flow velocity through the rotor (= ) 

/3 flapping angle 


<)c 7 _ 

daC 

doC oir ? c 

letting ^ - 

V / rJl 

ir - 

1 + \k sm 

/*l = 

^ / rJT 

= 

c 3 ^ 


) . + c i + s >- ^ 

a<r pod 


it comes : 


f>/; f ^ ^ 


dCz e + ^ 4 ^ C -p> / ) 


+ H tr 


cW 

Cz 


( y3 Cos V+o.'t,; + (Cj + S L a- } cf 


(-2; 


On other hand, the flapping equation may be written as : 


2 / J /I 
J* 2 


r -a 2 ( -ZU-± j 


K. 


■A 


r bdr 


C3; 


where fc> is the local chord, and I , the inertia of the blade 
about the hinge axis. 


It is assumed that : 


1) The induced velocity is uniform { = ete) 

2) The flapping equation can be approximated using the average 
value of the coefficient of /b (appearing in f r ) and ne- 
glecting the term in jj- fb Co s . 

Then Eq. 3 may be written as : 


d 2 A 


-2- J£— +/* 

SlT 0 


H 0 


) 
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2 

jxv 0 


A£z r 3 bdr 


(This equation defines the time constant ^ of the rotor) 
d M a - m 0 (VJ +- C?) u 


where R 


(V) *f[ 6+0- jLS? ( + fA.oiy) J rbdr (7) 


It is deduced that : 

1) /3 is obtained from M 0 by an ordinary linear transfer 
function. 

2) Mo results from the multiplication of the flap deflection by 
a function of V 7 . 

3) The aerodynamic forces may be expressed in the same way, 
i. e. , multiplication by a function of ¥ and then filtering 

by a linear filter. 


(F. *S^j — Fr 


Quantities associated with the rotor, such as flapping angle, 
vibratory forces, stresses and so on, which depend only upon 
the rotor dynamics, may be expressed as functions of the 
aerodynamic forces by means of ordinary linear transfer func- 
tions. Let X be such a quantity, and X ( co ) its Fourier 
transform, it comes : 

-R. 

X ( cu) - J G r (coj f r (co/ dr C?) 

- Such a model may be sketched as in figure 47. 

For steady state conditions, all functions can be expanded in 
Fourier series, such as : 

Jn j It 

x = 21 x n & 

tls ->> 

and a general form for x n is obtained : 


*g„ jf Sp <-p 



GIRAVIONS 
DO RAND 



Doc*. DH 2Q1X-JUES 


Page -22- 


This expression can be in turn transformed in a relation 
between the sine and cosine coefficients of the Fourier 
expansions of X and cP : 

X (V) s X 0 + x r c*s x's/n VS -t- x z cos 2 <V + X z sin ZVi- ... 
cP (V ) ■=, cP 0 + c*f cos VS -h cfj sin VS y- cos 2 *P -f c? 2 sin 2 Hs -+ 
constructing the vectors : 

x„ o J„ 

S X 0 , Ji 

x = x ’ x_ - ^'i cf = S, 


it comes finally : 

x =r TJ + x 0 c«; 

where T is a matrix independant of cf which consequently 
defines the relationships between cf , the input to the rotor, 
to X, considered as its output. 

When collective and cyclic terras ( cf 0 } <f 1 , c/ ; ) are constant, a 
similar relation can be written between the multicyclic com- 
ponents of X and cf : 

+ X Mo (12) 

This type of model was used in Ref. In addition, since the 
quantity in eq. 6 contained no higher harmonics than the 
third, it was shown that t m had some symmetry properties 
and could be written as : 


T = 
M 


2 

“ h 2 

g 23 

- h 23 

2 

g 2 

h 23 

g 23 

'32 

~ h 32 

g 3 

- h 3 

32 

g 32 

h 3 

g 3 


faasggag 
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2. 2 Results obtained with simulated data 

In Ref. 3 was described a graphical method for obtaining the 
coefficients g and h in T^. This was possible because of the 
particular choice of the flap deflection components. A general 
method has been developed later. It identifies all the coeffi- 
cients of T-yj without assuming any symmetry properties (see 
paragraph 3 below). The simulated runs, whose characteristics 
have been described previously in the introduction, have been 
processed by this method and the following matrix lias been 
found, relating the multicyclic part of the thrust coefficient 
to that of the flap deflection : 

C-r = T~ m v* Ct- o 

when cl 1 is expressed in degrees : 


47.2 

15.6 

- 0.1 

17. 1 

- 15. G 

47.3 

- 13.5 

- 8.0 

7.4 

- 15.2 

56.9 

12.5 

13. 8 

7.6 

- 13.4 

58. 1 


The symmetry properties found here are in quite good agree- 
ment with eq. 13, which is consistent with the fact that the 
simulation uses an uniform In flow field and that the variations 
of /3 are small. The residual lack of symmetry, however, 
results from the non-linear effects introduced by blade stall, 
reverse flow region and perturbations due to /$ . These de- 
viations are of course increasing with the rank of the harmo- 
nic which is considered. 

3. - A GENERAL METHOD FOR STUDYING THE ROTOR 
3. 1 Principle 

In steady state conditions, all variables of Eq. 2 and 3 may 
be expanded in Fourier series. The remarquable feature of 
these equations is that they are linear with respect to the 
variables /3, , but with periodic coefficients. This 

leads automatically to linear relations between the Fourier 
coefficients of these variables. 

Since the induced velocity can be computed in the vortex theory 
from the local aerodynamic forces, its Fourier coefficients are 
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also related linearly to those of fr . Eliminating fr between 
Eq. 2 and 3, will finally lead to linear relations between the 
coefficients of fb , or ^ , and those of . 

Finally, the Fourier coefficients of a quantity ^associated with 
the rotor can be linearly related to those of cf and Eq. 11 
remains valid : 

x « T$ + x Q 

But now, symmetry properties are no longer guarantee. 

Except for stall effects, this model should represent the present 
rotor with a good approximation. The problem now is to identi- 
fy the coefficients in T and x. from the measurements of 
X and cf obtained during the tests. 

Obviously, T depends upon p. , ~n. , £03 and ( on c^ s ; f 
and X and cf should correspond to constant values of these 
coefficients. 

3. 2 Linear least square multiple regression for multicyclic transfer 
matrix identification 

To proceed conveniently with this identification, it is desirable 
to write Eq. 11 in the more condensed form : 

X « T J 

where the last column of T is in fact the vector and the 
last component of cf is equal to 1 . 

Calling N,the number of Fourier coefficients of cf ( N = 2 x 
number of harmonics + 1), the vector cf has now M = N +1 
components and T a N x M matrix. 

Letting L be the number of sets of values ( X c5 ) correspon- 
ding to L test points, an M x L matrix fXl and an M x L 
matrix [cf] can be constructed by juxtaposition of the L columr 
vectors X and <S , and Eq. 7 may be written as : 

[X] = T [c?J (15) 

When L = M, T can be obtained by a simple inversion : 

t =. [x] [J] - 1 tie; 
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In such case, if a vector * is computed from a measured c? 
and the computed T given by Eq, 16, its value x is such that 

X - X cs o 

When L is greater than M, there is no solution to Eq. 8 
because in general, X and <S do not verify Eq. 14 exactly 
and Eq. 15 is a system of uncompatible equations. However, 
a value of T can be found which minimizes the quantity : 

<?=Z (X l _X ; ) T rX i -Xj= X CTci , ;.~X- l j' r (TcJ , ;.-X;J ( 17 ) 

(where the upper sc r ipt T denotes a transposition). 

Q is the quadratic error between the measured values X" and 
the computed {or "estimated’') values x . If x and cS verify 
Eq. 14 and no measurement error is made, Q is zero. Other- 
wise the value of Q is a measure of the correlation of X and 
(J with respect to the model. It can be seen indeed, that 
it would be approximately equivalent to choose M test points, to 
determine a value of with Eq. 16, and then apply Eq. 14 
to the L - M remaining test points. Letting : 

5 = £1 *; r X; 

l= r,i_ c u 

a correlation coefficient can be defined as : 


C 1 - 


(O/S) 


Since Q is a quadratic function of the "output" x of the model, 
its minimization leads to a least square estimate of T given 
by : 

T. IX) [ef] T ([JH«nV 

OnceT has been obtained by Eq. 19, the estimated X is com- 
puted from the measured tS by Eq. 14 and compared to the 
measured X for each test point i, in order to obtain : 

- the quadratic error 

a: = c xc - x-J T C x c ~ x c ) u o) 

- the relative error 

/ txjxi) < 21 ' 

- the averaged error is then defined as 

< e * = -t S )L 
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3.3 The problem of Ideal and optimal cams 

The main interest erf a multicyclic control of the rotor is the 
possibility of reducing blade stresses and hub vibrations by an 
appropriate choice of the multicyclic cam components. In the 
case of the jet flap rotor, the flap deflection can be separated 
into two parts : 

- the usual rotor controls (collective and cyclic pitch) which 
can be represented by a component vector, c/p . 

- the multicyclic control, represented by a vector containing 
the harmonics, d*? . 


Let’s consider for instance the case of the stress at some 
location on the blade. 

It is possible to partition Eq. 14 in the following way : 


<o = T J 



X 


* x x v v X 


X 


r 

V 


X X X 1 V X . . - X | X 

" 

Sp 



X 


vxx'xx x | x 


X 



' 




X 



• 


• * * 1 • . . t . 


X 

X — 


. 



• • * 1 • . 1 . 

X 

. 

1 



* T p j * 1 rn • 


4 , 



' 


• • * . • ' * * j * 


* 



X 


X X X X | x 


X 





^ > 


1 


So that the variable part of the stress may be written as : 

<Tm = Tp c? p + T m cf m + < 3 ^ C 23 ) 

It is therefore conceivable to act on(T m by changing cfln , and 
for given values of c? P and <&m a , find a value ofdVn which 
will minimize 6m in some way. 

In Ref. 3, the ideal cam was defined as the cam cancelling 
harmonics 2 and 3 of the stress. To be consistent with this 
definition, the ideal cam could be defined in general as the cam 
concelling the harmonics- 2, 3, ... NHARM of the stress. 
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In this case (N - 3) components of Cm have to be cancelled, 
and (N - 3) components of are available which can be 
therefore uniquely defined by solving Eq. 23. The corres- 
ponding multicyclic cam obtained in this case was denoted 
"A” type cam in Ref. 11. 

However, as far as the stresses are concerned, it is usually 
considered that the important reference parameter is the 
peak-to-peak value. Therefore it is more appropriate to deter- 
mine a multicyclic cam which would minimize the peak-to-peak 
value of the stress. Unfortunately, such an optimization problen 
cannot be solved analytically, but only on the computer, by 
iterative algorithms. 


As shown in Ref. 11, the optimal ’’A" cam does not necessarily 
reduce the peak-to-peak values because the decrease in har- 
monics 2, 3, ... magnitude may be associated with an increase 
in harmonic 1 magnitude which will result in an overall increa- 
se of the peak-to-peak level. Although not equivalent to a peak- 
to-peak minimization, the minimization of a quadratic function 
of all harmonics of the stress would give a closer, and there- 
fore better, solution. 

A type "B" ideal cam is defined as minimizing the quantity : 

^ oi — <5 m 

In this case indeed, G'm contains (N - 1) components and they 
are only (N - 3) in efm . Therefore there is in general no 
for which <T m = 0. 

However it is possible to find a value of which minimize 
Qm > and this value is given by : 

cf.-d. = - ( T m r 1 ( r p + <y mo ) t24) 

This shows that the ideal cam depends upon the rotor control 
. When several configurations are considered it is pos- 
sible to define an averaged value of . This average defines 
a cam which will have stress reducing effects in all cases, if 
the term T p c?p does not vary too much from one configuration 
to another. Such cam is called "optimal" and the "B" type is 
defined by : 

)' 1 T m ' r (T p < cfp >+G- m „) C 2S) 


where the symbol < > denotes an averaged value. 
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3.4 Fixed stick and fixed flight conditions 


Equations 24 and 25 define cams for fixed stick conditions 
since c?p is supposed to be unchanged. The stress reduction 
may be defined by comparing the stress peak-to-peak value 
when no multicyclic control is applied, with the case where 
the cam is introduced. 

An other possibility should be considered because, as it will 
be shown later (see paragraph B III), the forces acting on the 
rotor hub (L, X, Y) depend upon the collective and cyclic 
pitch, which is well known, but also upon the multicyclic com- 
ponents. 

Therefore the introduction of a cam will always modify the 
equilibrium of the rotor. In a real flight, the imposed con- 
ditions are lift, propulsion and lateral force, rather than the 
’stick” position (in fact the collective and cyclic controls). 

Therefore the control dp should be changed to maintain L, X 
and Y the same. 

Because L, X and Y can be also linearly related to the deflec- 
tion <T (as shown in B III), the stress can be expressed by an 
equation similar to Eq. 23 : 






where is a vector of components I/, X, Y. 

In the same way, the ideal cam for fixed flight conditions is 
obtained by : 

( ^LA ) ff c = - ( T m ) T fn ( T p G "'ol ( 27 ^ 


The method described in paragraph 3.2 readily applies for 
determining the matrices T*p , T ^ and the vector which 

can be casted into a single matrix T (like 1^, T m and <T m 
In T ) while a vector 3 is constructed with L, X, Y as the 
three first components, the others being the same as those of 
3 (i. e. multicyclic terms and 1.). 


It comes then : 


6 ~= T'cF 


/ 

and T' is obtained from measurements of 6 and 3 by equa- 
tion 19, which gives in this case : 

[s-] [j'][<J , 'cf' T ] -1 
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IMPROVEMENTS AND EXTENSION OF THE METHOD 


Matrix regularity 

Since the computation of T by Eq. 19 requires an Inversion, 
it is important that the matrix ([<$ ] [3 ] T ) be not only regular 

but also well-conditioned. Otherwise instability may appear 
in the coefficients of T . Ill-conditioning comes from a "linear 
closeness" between some vectors c? and its detection would in 
general require special algorithms. (See Ref. 5 for instance). 

In the present case, linear closeness occurs because the collec- 
tive value is always close to 24 degrees in all runs, and the 
last component of S is already equal to 1 in all runs. To 
avoid instability, 24 degrees is taken as a reference value 
for the collective, and it is the difference between the actual 
collective and this value which is entered as the first component 
of S . Also, before inverting the matrix, a quick check is 
made to eliminate pairs of linearly close vectors or null 
vectors. In this case some columns of T are missing since 
they cannot be determined from the data (this happens indeed 
with simulated data where cyclic and collective terms are 
strictly constant). 

Effect of the rotor shaft angle • 

Because of the limited number of runs available at constant ^ 
and o' a correlation analysis on more runs was not possible 
unless the effect of o^ s could be included in the model. 

This can be done simply by assuming* that the effect of o< 5 is 
independent of cf and can be represented by a vector <5^ such 
that Eq. 14 is replaced by : 


G = T3 


*3 ^5 


By adding to T an extra component representing 6^ and an extra con 
p one nt to ^ equal to ots , Eq. 14 can still be used and the same 
method remains valid for determining T by a least square 
regression. In this case however : 

M = N -» 2 

and, in the computation cf ideal and optimal cams, is to 
be replaced by : 


* This can be deduced from Eq. 2 since or q = + /2» 1 



QTOOA y- l YOU 





GIRAVIONS 
DO RAND 


► DII 2011 D E5 



4. 3 Effect of the advance ratio 

The variation of the rotor response with the advance ratio, 
at constant RPM, may be described by the evolution of the 
matrix T when p- varies. 

If T is known for a given value of jc , say , a first 
order expansion of T may be written as 

T < fM = T iy 0 ) ) Tj,. 

The matrix T"^ may be identified by the same method as T 
if the following transformation is made : 

= <5 - r CKo) 3 

3 ' = ( \A C? 


since 


& ( y ) - ~T ( n- ) <$ 


it comes : 


s' a Ta S ' 


( 


4. 4 Evaluation cf the stress reducing effect of the cams 

In the present tests, the mechanical cams have been tried with 
various phase and amplitude settings. As shown in Section C, I, 
stresses were reduced in some cases, increased in others. 

To evaluate properly the effect of ideal and optimal cams, 
stresses are first computed with no multicyclic cam and then 
compared to those obtained with ideal and optimal cams. If 
the ideal cam gives the maximum reduction, a slightly diffe- 
rent cam would obviously give a little less. To have a measure 
of the difference between the ideal cam and another one (like 
optimal or real), the following quantity is computed, "relative 
distance" between the cams in the Fourier space : 

cl - ^ ~ ^ C qL - <^rr> ) ( 30 ) 

where the vectors cf id and contain respectively the mul- 
ticyclic components of the ideal and the other cam. 
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5. - COMPUTER PROGRAM FOR T MATRIX IDENTIFICATION AND 
OPTIMAL CAM EVALUATION 

A FORTRAN IV program was written originally for the IBM 1130 
and then extended for the IBM 360/75. 

5. 1 Operations performed 

This program analyses a maximum of 30 test points given by 
the Fourier coefficients of the jet-flap deflection and the consi- 
dered stress, up to the 4th harmonic. It performs the following 
operations : 

a) Computation of the T matrix according to Eq. 19, 

b) Computation of the ideal cam for each test point, and of the 
optimal cam. Type "A" cams were used in a first version , 
of the program and are now replaced by type "B", 

c) For each test point 

cl Computation of the Fourier coefficients of the stress 
in two cases : 

from the measured S , 

. from the actual rotor control and the optimal cam. 

c2 Computation of the quadratic error Qi and the relative 
error (from Eq. 20 and 21), 

c3 Computation and comparison of the peak-to-peak values 
measured directly 

computed from the real flap deflection 
computed from the real efp with no cam 
computed from the real c$p and the optimal cam 
computed from the real cfp and the ideal cam 

c4 Computation of the relative distances to the ideal 
cam of the real and optimal cams, 

c5 Computation of the maximum flap deflections obtained 
with the real S P and the optimal cam, 

c6 Plot of the time histories (e. g. values of the stresses 
vs. Blade azimuth position) of the measured, computec 
and optimal stress. 

d) At the end of the program, the relative error < £ > and 
correlation coefficient C. 

5. 2 General organization of the program 
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5. 2. 1 Main program 

In the main the principal arrays are dimensionned and the 

G D 005 GEN. ' P-lotou A y .|v QN 







GIRAVIONS 

DORAND 



Dnr t DH 2011 D E5_ 


._-a2r 


following input cards are read : 


First card 
Second card 


Third card 
Next cards 


comments 

program parameters (number of runs, 
of harmonics, value of the collective 
reference) 

program options (variable , variable ja) 

contain the data concerning the runs 
to be analysed. For each run : 

. the Fourier coefficients of the stress 
(or any other quantities which need 
to be analysed) with eventually the 
value of . 

. the Fourier coefficients of the jet 
flap deflection with eventually the 
value of . 


When the jx option is used, an additional set of cards con- 
taining the matrix T* and the corresponding value of ja. 
should be read. 

5.2.2. Subroutine MUTMA (Multicyclic transfer Matrix Analysis) 


It constitutes the main core of the program and performs 
operations a, b, c, and d described previously. Options 
are available making possible : 

- to read in the matrix T (instead of computing it) 

- to compute the matrix T only 

- to read in an optimal cam 

- to perform the analysis without computing ideals and 
optimal cams. 

Several subroutines are used in conjunction with MUTMA 
performing special operations : 

- ahgfi 

Computes the amplitude and the position of the first maximu 
of each harmonic c£ a cam (given by its sine and cosine 
coefficients). 

- CAM 

Computes and stores the ideal and optimal cams. 
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- MATRA 

Computes the transposed of a matrix. 

- MINV 

Computes the inverse of a matrix (Internal computations 
are made in a double precision arithmetic). 

- PEAK 

Computes the values of a function given by its sine and 
cosine Fourier coefficients, at every five degrees, stores 
these values (can store 3 functions, e. g. measured, com- 
puted and optimal stress), computes the maximum, the 
minimum, the corresponding values of ^ , and then the 
peak-to-peak value. 

- PLOT 

Plot and print the values stored by peak 
A simultaneous plot of three functions is possible. 

- PROD 

Computes the product of two general matrices. 

- SMPI 

Check the regularity of a square symmetric matrix and 
proceed with a partial inversion if dependent or null vec- 


re found. For 

instance 

the 

matrices 

0 

A 


ab> 

d 

0 

0 

A 


a 

d/fc 

0 

0 



0 

0 

c 

0 

0 0 


0 

0 

0 

0 

A" 


0 

tya 

0 

0 


0 

0 

Vc 
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The listings of these programs are available at 
GIRAVIONS DORAND but will not be joined to the present 
report. However, an example of output will be presented 
in Appendix I. 

5. 3 Possibilities of the program 

It is worth noting that, due to the particular form cf Eq. It), 
each component of the vector X, no matter what it does repre- 
sent (i. e. whether it is actually a Fourier coefficient of a 
quantity or any other static value) can be separately related to 
the Fourier components of the flap deflection and °< s (i. e. 
each line of the matrix T is evaluated separately). 

Therefore this program can be used to correlate different 
quantities, independent of each other, to the flap deflection. 

A remarquable application of this possibility has been the 
correlation of the aerodynamic coefficient Clr , , Cy , etc. 

to the O and <*f s as described later in III, 6. 

Also, once T has been determined, it represents completely 
the rotor (for given ja, & and -TL ). 

Thus it is possible to predict the behavior of the rotor for 
arbitrary flight conditions or cam configuration. 

III. - APPLICATION TO THE TES . DATA 


1. - JET FLAP DEFLECTION FOURIER ANALYSIS 

In order to perform the multicyclic transfer matrix analysis, it 
is first necessary to know the actual flap deflection Fourier coef- 
fic ients. 

A total of thirty test points have been studied corresponding to the 
conditions : 


» o,4 _a = 250 rpm © o7 = 5* 

A first analysis was made at GIRAVIONS DORAND using a mechani- 
cal Fourier analyser, directly from the oscillograph records. This 
is not a very precise method, and later on, the harmonic analysis 
of the same 30 signals was made at NASA, AMES Research Center, 
using discrete measurements of the oscillograph records and a digi- 
tal computer process ing. 
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Values of the Fourier coefficients obtained by this method are given 
in Table VIII, up to the 4th harmonic. 

2. - ANALYSIS OF THE FLAP BENDING STRESS 

Signals corresponding to the flap bending stress at 0. 45 and 0. 7R 
have been first processed at GIRAVIONS DORAND from the original 
oscillograph records, in order to obtain 60 discrete numerical 
values per period. These sets of values vere then Fourier analysed 
on the digital computer at AMES. 

Fourier coefficients of the flap bending stress at 0.45 K are 
f iven in table DC Note that three runs were entirely analysed at 
GIRAVIONS DORAND. 

After encouraging results obtained with fewer runs (see ref. 6), the 
matrix analysis was performed for the 30 runs, first using the 
G. D. values for the flap deflection, then the NASA values. 

The error patterns obtained with G. D. and NASA values are dis- 
played in Fig. 48 were the errors £ are sorted by increasing 
values. They show that a good correlation is obtained in both cases 
for the 30 runs and that this correlation is increased with increa- 
sing precision in the measurements of J , indicating therefore that 
the model is correct but measurement errors are responsible for the 
scatter of the results. 


3. - 


To obtain a better idea of the statistical distribution of the error, 
the partition function has been plotted for both cases on Fig. 49. 

This function, F {£), is the ratio of the number of run whose 
error is less than £ , to the total number of runs. 

CORRECTION OF PHASE ERRORS 


Even with the NASA values, some runs still remained with a large 
error. For example run 16.08 exhibits a 20% error. If one looks 
at the waveforms of the computed and measured stress, as sketched 
in Fig. 50, it appears that the shape is almost correctly reproduced 
except for a phase shift, which in this case would be of about 10 
degrees. Because this error seemed of systematic nature, rather 
than the result of random errors on the computed or measured Fou- 
rier coefficients, the oscillograph records were examined carefully 
and it was soon discovered that the Ip pip, used to define the rotor 
period, was not regularly spaced but was under go ing random devia- 
tions about some average position. On the contrary, the period deduced 
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from the occurence of specific features of the signal itself was per- 
fectly regular. This is believed to be the result of some backlash 
occur ing in the slip ring drive. 

Since an error of l mm corresponds to 6°, and the undetermination 
on the lp pip is, in this run of one or two millimeters, it is not 
surprising that an error of 10° could have been made. 

In order to correct these errors, all the results concerning the 
30 runs have been systematically examined, and the phase shift 
estimated from the comparison between measured and computed wave 
forms. Then the Fourier coefficients of both <5 and have been 
modified consequently. 

These corrected values are given in Table X for J and Table XI 
for <5 o. 45 - The value of the correction^ is indicated in both tables. 

FINAL RESULTS FOR JET FLAP DEFLECTION AND FIAP 


BENDING STRESS 

With the values given in Tables X and XI, a matrix analysis has 
been performed. A very clear improvement of the correlation has 
been obtained as it appears on Fig. 48 and 49. 

It is interesting to compare the results obtained in different cases of 
study of the stress at 0.45 R, for fixed stick conditions, i. e. 
the input is only the Fourier coefficients of the iet flap deflection. 


Study Conditions 

cte 

G. D. values 

= variable 

G . D. values 
NASA values 
Phase cor- 
rected 




Number of 
test points 

Average 
error % 

Correlation 

% 

15 

6. 1 

81.6 

30 

9. 2 

85.5 

30 

8. 9 

86 

30 

6. 9 

89 


This shows that the correlation stays almost constant although the 
number of test points processed is doubled, which is an excellent 
proof of the statistical validity of the model. The effect of decreasing 
the measurement error is an increase of the correlation, which 
is also supporting this model. 
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4. 1 Fixed stick conditions 

The analysis of 30 test points from the values given in Tables X an< 
XI leads to a value of the T matrix elements given in Table 
XII. 

Boxes contain the coefficients of one harmonic upon an other 
one. In order to emphasize the role of the different multicy- 
clic components upon the harmonic content of the stress, each 
box may be reduced to one value expressing an average ratio 
of amplitudes. This gives in this case the following matrix. 
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Op 

f P 

2P 


4p 

Op 

1. 

0.2 

0. 15 

0. 13 

1. 

'p 

0. 33 

0. 71 

0. 08 

0.51 

0. 85 


0.41 

0.47 

1. 1 

0.34 

0.64 

3 P 

0.26 

1.1 

1. 2 

3. 9 

4.37 

Ap 

0.34 

1.37 

0. 72 

2. 1 

5.6 


where the coefficients have been normalized to • 

< o c ( = 0. 12 hbar/degree. 

This way of depicting the multicyclic effects is rather qualita- 
tive but gives a fair approximation of the influence of the iet 
flap deflection harmonics. It can be deduced for instance that 
1 degree of harmonic 4p in flap deflection will result into 
0. 85 x 0. 12 = G.lhbar amplitude of harmonic lp, 

0. 64 x 0. 12 = 0.072 hbar amplitude of harmonic 2 p and so on. 

These coefficients show also the remarquable increase in effi- 
ciency as higher harmonics are concerned, which could be the 
result of the rotor blades being near a first mode resonance 
for frequencies corresponding to harmonic 4p. 

This is emphasized by plotting the harmonic amplitude of the 
stress as a function of the rank of the harmonic of S .In 
general & n - P <^n) would represent on ordinary linear trans- 
fer function. Due to non linearities <j~ n => f (cf p ; is not zero and 
may be interpreted as cross transfer functions. These curves 
are displayed in Fig. 51. 
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In fig. 52 are represented the stress harmonic sensitivity to 
rotor shaft angle and flap deflection harmonics. The nu- 
merical results of the analysis are given in Table XIII. As 
mentioned previously the correlation obtained is fairly good, 
the average error being of 6. 9?< . 

The different ideal "B" cams and the optimal cam correponding 
to those 30 runs are sketched on Fig. 53. It is remarquable 
to note that the phase dispersion is small for the ideal cams, 
but amplitude may be quite different. The dispersion in the 
ideal cams comes from the different values of the control Jp 
and angle c/ s according to Eq. 24. 

If T p <S p is not varying too much compared to i^m 0 + )» 

then the dispersion is small. 

The optimal cam found in this case is the following : 

J ' l = 3. 9 Cos 2 Y - 1. 1 Sin 2 ^ - 3.1 Cos 3^- 0. 1 Sin 3f - 

2. Cos 4 ^ - 0.9 Sin 4 

In Fig. 54 are plotted the stress reduction obtained with the 
optimal "B" cam and the real cam (normalized to the ideal 
"B" cam reduction), versus the ’'relative distance" defined by 

Eq. 30. When this distance is small, the reduction obtained 

becomes comparable to the ideal reduction. However, since 
the ideal cam does not exactly minimize the peak-to-peak values 
in some instances the optimal cam could give better results, 
than expected. Nevertheless the general pattern indicates that 
the quadratic criterion is a good approximation for rapid in- 
vestigations. If a type "C" cam is defined as the cam mini- 
mizing the peak-to-peak, type "B" could be used as a starting 
point in an optimization algorithm for determining the "C" cam. 

4. 2 Fixed flight conditions 


The results concerning the fixed flight conditions are of course 
of greater interest since it corresponds more closely to the 
real helicopter case. 

To run this analysis, the values of ^ and c?, have been 

replaced respectively by those of /g“ , C XR /$- and Cy /fT 
extracted from Table II. 
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(In fact, these three coefficients have been multiplied a 
factor I0 3 in order to have reasonable coefficients in the 
matrix T). 


The corresponding matrix T found in this case is given in 
Table XIV. It differs from the matrix obtained in fixed stick 
conditions because of the influence of the multicyclic flap 
deflection upon the values of C L , C x and Cy. As will be 
seen in paragraph 6, 




c lr / ( r 
|c*r /g- 

Cy/6- 


= r p cf p + F m + F f 


Eq. 26 gives then 

= (Tp + Ct; F, 

comparing to Eq. 23 it 


TY\ 


m. 


. + V 


,) 


comes 


r P 


Vo 


However it seems that the term 7^'^ is not very important 
as attested by the similarity between T m and (e. g. see 
Tables XII and XIV). 

Results concerning the analysis of the 30 runs are given 
in Table XV and deserve some comments. 


Although the error is 2% higher than in the fixed stick case 
(due, it seems, to some errors in the measurements of 
Cls » Cx * and Cy ), the match between the model and the 
data is still good. 

The reductions obtained with the ideal cam reach values 
as high as 66" and is never less than 40". The averaged 
reduction is 49%. 


With the optimal cam, an average reduction of the stress of 
36% is found. 


These results seem to agree very well with the reductions 
observed experimentally. 

Ideal and optimal cams are sketched in Fig. 55, exhibiting 
the same concentration of the phases. 

On Fig. 56, the optimal cam is compared with the experi- 
mental cam IV. 
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Except for the 2p harmonic, there is a remarquable similarity 
between these two cams. 

The optimal cam for fixed flight conditions is given analyti- 
cally by {fig. 62) : 

3.4 Cos 2 V + 1.4 Sin 2^ - 5.7 Cos 3H 7 - 4.5 Sin 3 V 
- 7.5 Cos 4^ - 4.8 Sin 4^. 

It is worth noting that the harmonics amplitude though greater 
than in the fixed stick conditions, are still reasonable and 
could be actually obtained on the jet flap. 

5. - ANALYSIS OF THE VIBRATIONS 

5. 1 Introduction 


The method used successfully for the analysis of the flap 
bending stresses can be equally well applied to the problem 
of the vibratory forces on the rotor hub. 

Reduction of these vibrations are indeed a major problem 
in helicopter technology and the multicyclic control of the 
rotor appears to be a very powerfull mean of controlling 
these vibrations, as already seen in Section B, I. 

The forces and moments acting on the rotor hub were measu- 
red by a 6 - component balance. The forces on the three 
vertical arms of the balance correspond to the signals recor- 
ded by the galvanometers 5.1 (left aft arm), 5.2 (right aft 
arm), 5.3 (forward arm). 

The sum of these signals gives the vertical force acting on 
the hub. The variations of this force define the vertical vibra- 
tions induced by the rotor, which will be transmitted to the 
fuselage in case of a real helicopter. 

5. 2 Processing of the test data 

The Fourier coefficients of the vibratory forces have been 
obtained from the oscillograph records using a mechanical 
Fourier analyser. The value of the coefficient is read on 
the instrument in millimeters, with an error of j_ 0. 05 mm. 
This corresponds to an error of 2 daN in the recordings 
considered here. In fact, because of the difficulty of following 
exactly the waveform, the error is about 6 daN. 
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The coefficients for the vibratory force have been finally 
obtained by summing up of the signals 5.1, 5.2 and 5.3. 

Because of the amount of time required, the analysis has been 
limited : 

- to the even harmonics (0, P2 and P4), since the odd ones 
should normally cancel out in a two-bladed rotor, if non- 
symmetrical effects are not present. 

- to the minimum number of runs needed to obtain a matrix T. 

For this second condition, 11 runs are necessary when the 
option is used and 4 harmonics considered. 

Using the dependency analysis described in Ref. 10,11 vectors 
^ have been chosen among the 30 considered in the analysis 
of the stresses in order to obtain a basis as orthogonal as pos- 
sible. A twelfth run was added to these to evaluate the corre- 
lation. 

In Table XVI are given the Fourier coefficients, for these 
12 runs, of the signals 5.1, 5.2, 5.3 and the resultant vibra- 
tory force F z . 

Analysis of the vibratory force 

When the multicyclic transfer matrix analysis is applied to the 
vibratory forces, a particular problem appears due to the fact 
that harmonics lp and 3p are absent. The matrix T has thus 
four rows of zeros as can be seen in table XVII (fixed stick 
conditions, f. s. c. ) and table XVIII (fixed flight conditions, f. f. c) 

In this case the matrix (T^T^ ) of Eq. (24) is irregular and 
the program normally fails when computing the ideal cams. The 
physical reason for this is that there are four quantities to 
define the vibrations (p 2 ,fg , f 4 , f 4 ' ), but six variables to adjust 
in the flap deflection (c^, t£, in order to minimize 

the vibration. 

The problem is therefore undeterminated. It is possible, however, 
to solve this by deciding not to consider harmonic 3p in the 
multicyclic cam. In this case, the vibration may be exactly 
cancelled out since its four Fourier coefficients are linearly 
related to the four Fourier coefficients of the cam. - 
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A slight modification of the program "CAM 1 ’ is necessary to 
compute the ideal cams. 

Its consists in setting to zero, the columns of corresponding 
to harmonic 3p. In this case the subroutine SMPI will reco- 
gnize corresponding zeros in (T m T 7Vy|) and proceed with the 
inversion of the reduced matrix without any problem. 

Ideal and optimal cams are plotted in figure 57 for fixed stick 
conditions and figure 58 for fixed flight conditions. It should 
be noted that the dispersion among the ideal cams is of the same 
order as for the stress cams. 

In tables XIX and XX are given the results of the analysis for ’ 
fixed stick and fixed flight conditions. The correlation is very 
good. 

The vibration reduction by the ideal cam is always 100°; but the 
optimal cam gives an averaged redaction of 40 7 (for f. s. c. ) 
and 29°.( (for f. f. c. ) which shows that the same type of cam 
could be used in all flight conditions (with h = 0.4). 

The optimal cams for vibrations are given by • 

f, s. c. 

24.1 Cos 2 H 7 - 0.6 Sin 2f + 7.2 Cos 4 'V + 11.9Sin4‘f r 

f.f.c. 

c£p*= -5.2 Cos 2H'- 8.2 Sin 2^ + 2.5 Cos 4^4 1.1 Sin 4^ 

It is interesting to note the difference with the cam obtained 
for reducing the stresses. Also these cams maybe compared 
to the optimal cam obtained in the simulations analysed in 
paragraph II 1 and II 2. 2. These simulations give the coeffi- 
cients cf c-r , which is closely related to the vibrations. The 
The cam obtained in this case (for fixed stick conditions) is : 

<S op t= -13.2 Cos 2<V-2. 8 Sin 24 ' + 0.3 Cos 3^ + 2. Sin 3^ 

and reduction obtained was 64 ?f. 


Harmonic 3p was considered here because non- symmetric effects 
were considered. In figure 59 are compared the simulation 
and real test optimal cams. 
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Finally, in figure 60 and 61 are represented the transfer 
functions and harmonic sensitivity as was done previously 
for the stresses, showing the influence of the flap bending 
mode resonance. 

ANALYSIS OF THE AERODYNAMIC COEFFICIENTS 


6. 1 Matrix T 

It is possible to apply the matrix analysis to correlate some 
aerodynamic coefficients to the flap deflection and rotor shaft 
angle oL . 

The same program is used, but the aerodynamic coefficients arc 
introduced instead of the Fourier coefficients of the stresses 
(or vibrations). 

In this case each row of the matrix T gives the derivatives 
of the corresponding coefficient with respect to cP and . 

The coefficients which have been anafysed are : 

C L r / ^ rotor lift coefficient 

/ G rotor propulsive force coefficient 
Cy / <o lateral force coefficient 
cm / <3 pitching moment coefficient 
Cr r / roll ing m om ent c oeff ic i ent 

As before, the same 30 runs have been studied. 

The corresponding matrix T is given in table XXI. 

The influence of mult icy cl ic inputs is quite remarkable espe- 
cially concerning the 3p and 4p components. 

Contribution to lift of these terms is important. For instance, 
the optimal cam (for fixed flight condition) for reducing the 
stresses, contributes to about 20% of the lift. 

6. 2 Flight control matrices 

Besides the role of the multicyclic components which is still 
difficult to explain, it is interesting to know the conventional 
control effects on the rotor forces. 
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The constant term in the Cl. r/g* coefficients comes from 
the blade pitch angle & q ? which provides a lift even 
for no flap deflection. By construction of the blade, 
the lift is zero when 6 ^ ^ - 4°. 

It comes therefore : 

f0 3 C L * (G = 3.3 J c + °- 2 + 12 - * 6A 

For the propulsion : 

TC 3 Cxr i <3 *0Z<S o +0.2c -it .3 

Here the value - 11.3 represents the drag of the rotor. 

It is interesting to note that for ~ -10°, this drag is 
just balanced by the rotor tilt. 

For the lateral forces : 

I0 5 Cy /CT = 0.3 <$ 0 - 0,3 <?, + o.l e? 1 / - 0.9 4- 0.4^5 

The value - 0. 0 indicates that the rotor tends normally 
to tilt to the left (when facing the wind). This increase mar- 
kedly when the rotor shaft is tilted forward. 

From the three preceeding equations one can deduce 
the controls to be applied to the rotor in order to obtain 
given forces : 


r S 0 - 2,67 (f0 3 C LR /G -107.6) + 0.248 (to 3 Cx*/<S 4-11.3 ) + 0,16 (T0 3 Cy /<5 + o,9).2W< t 

. cf 1 - 0.583 C 10 3 C t ft /<S - <07. 6) _ 1.48 C to 3 C xR /<S +ii .3 ) -4.9 C™ 3 Cy /^ + ^)_4 5<* 5 

= 0. 46 O0 3 Ct. R /(5- -To7.8) _ 3.27 ( <0 3 C„«/<? + n .3 ) _ 2.39 C f0 3 C y /S'- o.9)- 


(It is to be recalled that these coefficients are only valid 
for f*. - 0.4, e o 7 = 5° and jT = 250 r.p. m.). 

6.3 Analysis of Performance characteristics 

The domain of the tests, presented in tables 1. 1 through 1.23, 
comprises 29 runs and 334 test points performed at wind 
tunnel speeds ranging from 0 to 123 knots, highest speed 
attained (see figure 63). These basic data has been processed 
by standard computer techniques, tables II. 1 through 11.24. 
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The results obtained have been presented graphically, see 
figures 1 through 20. 

The tests have been focused on comparatively high tip- speed 
ratios attaining 0. 7. However, neither maximum nor optimum 
capabilities of the jet-flap rotor could be demonstrated during 
this series of tests because of lack of power of the compressed 
air supply generator. The comparison with the first series of 
tests shows that the performance remains unchanged since the 
modifications introduced in 1965, fig. la. The test domain has 
been enlarged to very high blade loadings, the highest 
recorded in hover configuration reached 3.26 for a ratio of 
C£ - R /c.j R of 4.25. The highest Ct -» ? /Cj , R ratio, which is com- 
parable with the *-/ D ratio, is 6. 5 at Cl equal to 1.07. At 
advance ratio 0.25 we have Cl -R /cj R = 6.1 at = 1.05. With 

increased tip-speed ratios these characteristics are somewhat 
lower, at V/SLR - 0. 6, C <~R /Cj R = 5.1 at = 0.97. At the 

same time high propulsion forces are maintained at forward 
flight. At advance ratio of 0. 6 the ratio L is still 12.4%. 
As we mentioned these characteristics are neither the highest 
possible nor the optimal. A better power adaptation would permi 
a much more effective demonstration of jet-flap rotor capabili- 
ties. 

As seen in the figures 1 to 20 a high degree of scatter existed 
during these tests. In fact difficulties were encountered to 
maintain all but one parameter constant during the tests. In 
particular, secondary effects such as multicyclic effects have 
been present in most tests. It would be useful to employ a 
more advanced technique of analysis to obtain more consistent 
results. Such an approach has been outlined above but has not 
been per sued much further because of limited time allotted to 
the present analysis. 
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CONCLUSIONS 


1. DATA PROCESSING 

The 1971 tests of the Dorand H-2011 rotor have been analyzed to answer the 
following two main themes of research : 

a) Effectiveness of multicyclic effects on the reduction of stresses and vibrations 

at high advance ratios. 

b) Influence of minor modifications which the rotor has undergone since the pre- 
vious tests, in 1965, on the general performance characteristics. 

On the outset of the data processing work it became obvious that a special 
analytical procedure was necessary to obtain a clear picture of the multicyclic 
phenomena we encountered in the experimentation domain covering some 300 test 
points and ^u values ranging from 0.25 to 0.70. Such a procedure lias been establi- 
shed and has permitted to deal with the total experimental data output in what can 1 
considered a simple and global manner. The basic idea behind the procedure used 
consists in relating the output values, stresses or vibratory forces, to input signals 
multicyclic jet deflections, by a transfer matrix. Analytical and experimental work 
has shown that such transfer matrices exist and are practically invariant for given 
advance ratio. This procedure, it is believed, opens a very wide field of investi- 
gation to most rotary wing problems concerning stresses and vibrations. 

The existence of a unique matrix representing the transfer of control inputs 
to measured stress and vibratory forces has permitted to optimize the multicyclic 
components to obtain maximum stress or vibration reductions. The test- have 
already shown reductions reaching 50% of stress or vibration levels. These reduc- 
tions have been obtained by a sweeping method, A more systematic approach during 
the first series of tests has not been possible for the following reasons : 

1) The data obtained could not be analysed instantaneously. 

2) No simple method existed to guide the choice of multicyclic laws to 
obtain better reductions of stresses and vibrations. 
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3) The correlation between the control inputs and multicyclic effects of the 
jet flap was comparatively poor, because of mechanical wear and back- 
lash in control linkages progressing rapidly with number of runs. 

2. VIBRATORY STRESS ANALYSIS 

The stresses measured at 0.45 and 0. 70R have been recorded and Fourier 
analyzed for 30 test points at ^u 0.4. Because of near resonance conditions the 
results show a strong effectiveness of 3 P and 4 P signals e. g. a 4 P jet flap 
input results in a 4 P stress 5. 1 times higher that 2 P stress resulting from 
a 2 P jet flap input. The interharmonic couplings appear very strong. The ana- 
lysis has shown that stress harmonics higher than 5P can be neglected. The 
sensitivity to the jet flap deflections varied from 0. 10 hbar per degree of jet 
flap deflection for 1 P and 2 P to 0.48 for 3 P and 4 P control inputs. Similar 
trends are shown in the study of the bending stress sensitivity to the rotor shaft 
angle. The correlation obtained is fairly good, the average error being 6. 9%. 

Though the tests have been done in "fixed stick conditions", the "fixed flight 
conditions" which are of greater interest as they correspond more closely to the 
real helicopter case, have been investigated. The stress reductions obtained with 
"ideal cam" reach values as high as 66% and never less than 40%. In both cases 
the match between the matrix model and the data is good, though the error is 2% 
higher for the fixed flight conditions. 

The computational difficulties to use the peak-to-peak level as an optimization 
criterium have necessitated the introduction of multicyclic laws ("ideal" cams and 
"optimal" cams) reducing in some quadric sense the harmonic content, but not 
necessarily the general peak-to-peak value. It was remarkable to find that the 
ideal and optimal cams exhibit the same phase characteristics. Moreover, when 
compared with the optimal cam, the experimental cam N° IV shows a noticeable 
similarity except for the 2 P harmonic. 

The optimal cam for fixed flight conditions is given analytically by : 

Sopt = 3.4 cos 2 V - 5.7 cos 3 V' - 7.5 cos 4 V" 

+ 1.4 sin 2 f - 4.5 sin 3 Y - 4. 8 sin 4 
Sopt being expressed in degrees of jet flap deflection. 
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3. VIBRATORY FORCES 

The vibratory forces on the rotor have been measured by load cells situated 
beneath the hub and Fourier analyzed. Only vertical forces of 0, 2P and 4P com- 
ponents have been taken into account. The same matrix method has been used for 
analysis of the bending stresses with the modification cancelling the 3P 8 com- 
ponents. However, to eliminate the undetermination stemming from the fact that 
the flap deflection was defined by six variables and the vibrations by four only, 
the analysis used 12 test points to evaluate the transfer matrix. 

The highest vibration reduction noted experimentally was 4 8% (run 25.2 and 

25. 5). Analytically defined optimal cam laws differ for fixed stick and fixed flight 
conditions. In general the multicyclic effects which are optimal for stress reduc- 
tions are not the same as for vibration reductions. 

The limited scope of work has not permitted to obtain a clear understanding 
of phenomena which govern the optimization of multicyclic laws. It is believed 
that a considerable research effort is necessary to define the cause - effect links 
and the sensitivities to various rotor parameters. It is also suggested to improve 
the investigation techniques and replace mechanical cams by electro- hydraulic 
devices which are more flexible in obtaining desired multicyclic laws. 

At ^u = 0.4 and for fixed stick conditions the optimal cam for vibration 
reductions gives an averaged reduction of 40%. The multicyclic law of the cam 
expressed in degrees of jet flap deflection is as follows : 

S = 24. 1 cos 2 y - 0. 6 sin 2 V + 7.2 cos 4 V + 11.9 sin 4 
opt ' 

Introduction of non- symmetric effects seems to improve the results. When 
based on dynamically variable C T coefficient, which is closely related to vibrations 
on the hub, the reduction obtained was 64%. 

4. PERFORMANCE AND AERODYNAMIC ASPECTS OF THE ROTOR 
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As already shown during the whirl tests outside the tunnel the modifications 
introduced since the 1965 tests had little influence on the performance of the rotor. 
The C TT1 /cr against C TTJ /<r curve remained unchanged. The only noticeable 

LiK d JK 
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change that the modification (which consisted in the increase of the nozzle height 
from 4. 2 mm to 5. 3 mm and in the increase of the jet flap deflections from 
± 25° to ± 50°) have brought, was the shift of the zero lift chord, -4.0 degree 
nose down. 

Most of the tests have been performed below the optimal power conditions 
corresponding to the pressure ratio of 4. 0. In fact due to the power supply 
deficiencies the maximum pressure ratio obtained during the test did not exceed 
3.6. The rotor nevertheless produced very high lift and propulsive forces. The 
highest mean lift coefficient, t attained on the blades was 3.2. At advance 
ratio of 0.6 the propulsion force was still 10% of the lift for o(, s - 12°. The 
mean lift coefficient was 1.3. These figures are typical to the jet flap rotor. 

An attempt was made to correlate the aerodynamic coefficients to the flap 
deflections and the rotor shaft angle. The influence of the multicyelic inputs is 
quite remarkable. For instance the optimal multicyelic law for stress reduction 
contributes to about 20% of the total lift. 

5. FINAL REMARKS 

The main conclusions of the work done can be stated as follows : 

1) The multicyelic effects bear a very strong correlation to the stress 
levels and vibratory forces. 

2) The reduction of stresses and vibrations obtained by non optimized test 
procedures has reached 50% at advance ratio 0.4 (Maximum stress and 
vibration reductions did not, however, occur at the same multicyelic 
settings) . 

3) The performance characteristics of the rotor do not seem to be essen- 
tially modified by the multieycUe inputs, their contribution may, however, 
amount to as high as 20% of the total lift (This effect seems to be due 

to compensation of cyclic inputs by multicyelic inputs in flight fixed 
condition and has been noted in the domain far from the stall limitations). 
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4) New method of analysis of multicyclic effects has permitted optimization 
of control inputs. This method is based on transfer matrices, which are 
practically invariant with blade angle setting, shaft angle, and basic 
control settings for a given advance ratio. The transfer matrices indicate 
that, at advance ratio 0.4, for example, stress and vibration reductions 
of 80% are feasible. 

5) The modifications introduced to the jet flaps since 1965 tests did not 
influence the performance of the rotor to a noticeable degree. 

The limited scope of the test analysis did not permit to draw definite conclu- 
sions concerning the maximum attainable efficiency of the jet flap rotor. However, 
farther improvements reaching 80% reduction of stresses and vibration are theo- 
retically feasible. 

The matrix analysis used to correlate dynamic outputs to multicyclic inputs 
has been of great help in gaining an initial understanding of the complex relation- 
ships between the aeroelastic and aerodynamic phenomena. The work done has 
lead to the discovery of the existence of a transfer matrix that is invariable 
under many rotor configuration variations and describes in a most global and 
simple way the behavior of the rotor in forward flight. It is believed that this 
new approach to the analysis of the rotary wing phenomena can make a decisive 
step in solving the inextricable complexity of the mathematical modeling of the 
helicopter rotors working conditions in most flight configurations. 
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DE SCRIP TION OF THE COMPUTER PROGRAM OUTPUTS 


In this appendix is given an example of the results of the analysis by transfer 
matrix program. 

This example lias been taken from the analysis of the flap-bending stresses 
at 0.45 R. Thirty runs have been analysed ; the multicyclic cams are computed 
for fixed stick conditions. 

The results corresponding to one test point are given on two sheets of machine 
paper (tables XXII and XXIII). 

1. FIRST SHEET (see table XXII) 

1.1 Title and test point number : in our example stresses for the rim 
number 14 - 13. 

1.2 <^ s value (here = - 10 degrees) 

1.3 Four columns containing Fourier coefficients (OP, IP cosine, IP sine, 

2P cosine, etc. ) of the optimal, computed and measured stresses and of 
the difference between computed and measured stresses. 

Units for these values are given in the table (here hectobars). 

1.4 Then is given, the value of the correlation, the quadratic error and the 
relative error (5.1%). 

1.5 The following table indicates, for the computed, measured and optimal 
stresses, the value and the azimuthal location of the maximum and 
minimum cf these functions. It gives also the peak-to-peak value. 

1.6 After that, for four multicyclic inputs (i. e. no cam, real cam, optimal 
cam and ideal cam) are computed the following quantities : 

a) - The relative distance to the ideal cam which by definition equal 

0. , for the ideal cam and 1. for the "no cam" case. 

b) - The peak-to-peak stress reduction (or increase when negative) with 

respect to the peak-to-peak value when no cam is used. 

c) - The corresponding peak-to-peak value. 

1.7 Then is given the maximum and minimum of the jet flap deflection for 
the optimal cam and their blade azimuthal locations. 
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2. SECOND SHEET (see table XXIII) 

On this sheet are plotted the time histories of the optimal stress ( . ), com- 
puted stress ( * ) and measured stress (+). 

The stresses values for these points are written on the right part of the 
sheet. 

It can be noted that a priority has been established between the symbols , 

+ and . on the plot, 

if *■ and 4 (or . ) are at the same place, ★ is printed, 
if + and . are at the same plade, + is printed. 

Nevertheless, when this happens, there is no ambiguity about the location of 
these points because the three corresponding values are known. 

Similar outputs are obtained for each test point analysed, but we must remark 
that the program gives first the matrix T, the Fourier components of the 
ideal cams and optimal cam. 
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